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Recent developments in thermal leptogenesis are reviewed. Neutrino mixing data favor a 
simple picture where the matter-anti matter asymmetry is generated by the decays of the 
heavy RH neutrinos mildly close to thermal equilibrium and, remarkably, in the full non 
relativistic regime. This results into predictions of the final baryon asymmetry not depending 
on the initial conditions and with minimized theoretical uncertainties. After a short outline of a 
geometrical derivation of the CP asymmetry bound, we derive analytic bounds on the lightest 
RH neutrino mass and on the absolute neutrino mass scale. Neutrino masses larger than 
0.1 eV are not compatible with the minimal leptogenesis scenario. We discuss how the results 
get just slightly modified within the minimal supersymmetric standard model. In particular 
a conservative lower bound on the reheating temperature, Tr > 10 9 GeV, is obtained in the 
relevant effective neutrino mass range rh\ > 3 x 10 -3 eV. We also comment on the existence of 
a 'too-short-blanket problem' in connection with the possibility of evading the neutrino mass 
upper bound. 



"Compendium of'-' and'-' mostly based on'"' with some new results in 3.4, 3.8 and 3.9. 



1 Introduction 

Cosmic rays and CMBR observations indicate that our observable Universe is baryon asym- 
metric^. Moreover the observation of the acoustic peaks in the power spectrum of CMBR 
"■", combined with large scale structures observations^, provide a precise and robust mea- 
surement of such an asymmetry that can be expressed in terms of the baryon to photon 
number ratio at the recombination time, 

rg MB = (6.3 ± 0.3) x 1(T 10 , (1) 

in very good agreement with the latest determination from (NACRE updated) Standard 
BBN and primordial Deuterium measurements that give^ 

r? |^^ = (6.1±0.5)xl0- 10 . (2) 

At the same time there is a growing evidence that an inflationary stage occurred during the 
early Universe. In this case this would have diluted any pre-existing initial asymmetry to 
a level many orders of magnitude below the measured value, thus requiring an explanation 
of the observed baryon asymmetry in terms of a dynamical generation, the aim of a model 
of baryogenesis that necessitates the accomplishment of the three famous Sakharov's con- 
ditions: C and CP violation, B violation and departure from thermal equilibrium. Within 
the Standard Model all three conditions are fulfilled, yet the observed value is too large to 
be explained and therefore a successful model of baryogenesis requires some new physics 
ingredient. A host of models have been proposed since the first Sakharov idea^. Some 
examples of typologies of baryogenesis models are: Planck scale baryogenesis, baryogene- 
sis from phase transitions, Affleck-Dine models, baryogenesis from black holes evaporation, 
models of spontaneous bary ogen esis^. __„_ 

Even though leptogenesis^* and GUT baryogenesi s * * exhibit, from a particle physics 
point of view, substantial differences, they can be jointly regarded as two different examples 
belonging to the oldest class of models of baryogenesis from heavy particle decays. Such a 
classification privileges the thermodynamical aspect enlightening general properties that do 
not depend on the specific particle physics framework. We will thus discuss the kinetic 
theory of heavy particle decays in the first part, while in the second part we will see how 
leptogenesis is a specific remarkable example in which the new physics ingredient is provided 
by the seesaw mechanism and such that the observed baryon asymmetry is nicely related 
to neutrino mixing data. 

2 Baryogenesis from heavy particle decays 

2.1 Out- of- equilibrium decays 

Let us consider a self-conjugate heavy (Mj S> Mew) particle X whose decays are CP 
asymmetric, in such a way that the decaying rate into particles, T, is in general different 
from the decaying rate into anti-particles, f, and such that the single decay process into 



particles (anti-particles) violate B — L by a quantity Ab~l (—Ab-l)- The CP asymmetry 
parameter is then conveniently defined as 

^¥Tf (3) 

For a joint discussion of baryogenesis (Ab-l > 0) and leptogenesis (A#_£ < 0) models, 
it is useful to introduce the quantity e = Ab~l £• The total decay rate Trj = T + T is the 
product of the total decay width, Fg st , times the averaged dilation factor (1/7) 

rD=r rest/i\ (4) 

Sphaleron processes, while inter-converting B and L separately, leave B — L unchanged'^' 
and for this reason the kinetic equations get much simpler if the B — L evolution is tracked 
instead of the separate B or L evolution. Moreover it is convenient to use, as an independent 
variable, the quantity z = Mx/T and to introduce the decay factor D = Tb>/(H z). Another 
useful choice is to track the number of X particles, Nx, and the amount of the asymmetry, 
Nb-l, in a portion of comoving volume R 3 normalized in such a way to contain, averagely 
in ultra-relativistic thermal equilibrium, just one X particle (i.e. Nj?(z <1) = 1). 

The simplest case is when the X life-time, r = l/r^ st , is much longer than the age 
of the Universe, tjj = (2i7) _1 , at z = 1, when the X particles become non relativistic. 
In this way decays will occur when the temperature is much below the X mass and the 
X-production from inverse decays, or other possible processes, is Boltzmann suppressed. 
In this situation decays are the only relevant processes and the kinetic equations for the 
X-abundance and the B — L asymmetry are particularly simple to be written, 



dNx 

dz 
dN B - L _- dN x 

dz dz 



D{z)N x {z) (5) 

(6) 



and solved, 



Nb-l(z) = N B _ L + e[N x -N x (z)] (7) 

N x {z) = N x e~ r ^ dz ' D{z>) . (8) 

The solutions can be fully described just in terms of the decay parameter 

■prcst 

K = j^, (9) 

in terms of which D = K z (I/7). The dilation factor, averaged on the Boltzmann statistics, 

is simply approximated by'"' 

1\ z , . 

(10) 
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Figure 1: out of equilibrium decays. 
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a useful simple expression that makes possible to solve analytically the integral in the Eq. 
(jSJ), yielding the result 



N x (z)~N x e 



-K [4_1|* 



+ (f) 2 ln(l+A,) 



In particular the final B — L asymmetry is given by 

N B -L = Nh_ L + eN x . 



(11) 



(12) 



The baryon to photon ratio at recombination can then be obtained dividing by the number 
of photons at recombination N* ec (about thirty times the number of photons at the onset of 
X decays) and taking into account that sphalerons will convert only a fraction a sp h — 1/3 
of the B — L asymmetry into a baryon asymmetry. In this way one can write: 



IN* 



rjB ^ 



B-L 



3 iV^ 0C 



(13) 



It is useful to introduce the efficiency factor defined as the ratio of the asymmetry produced 
from the X decays, excluding the contribution from a possible initial quantity, to the CP 



asymmetry, i.e. 

N B _ L (z)\ Ni Q 
k{z) = 1 . (14) 

In the case of out of equilibrium decays one has k(z) = N x — Nx{z) and the Eq. Q can 
be re-casted as 

N b „ l (z)=N b _ l + ek{z). (15) 

This definition is such that the final efficiency factor, K{ = k(oo) = N x , is equal to unity 
in the case of an initial thermal abundance with zi <C 1. In Fig. 1 we show two examples 
of out of equilibrium decays, for K = 10~ 2 and K = 10~ 4 , assuming an initial thermal X 
abundance (N x = 1) and zero initial asymmetry (N B _ L = 0). The numerical results are 
compared with the analytic expression (cf. pijl). 

The out-of-equilibrium picture is an efficient way to produce an asymmetry from decays. 
However it relies on the possibility that an initial X abundance was thermalized by some 
unspecified mechanism at T > Mx and that one can neglect a possible N l B _ L generated 
during or after inflation and before the onset of X decays. Therefore, it is evident that this 
picture is plagued by a strong sensitivity to the initial conditions and hence it requires to 
be complemented with a model able to specify them, for example a detailed description of 
the inflationary stage. 

2.2 Inverse decays 

The out-of-equilibrium picture is strictly valid only in the limit K — ► 0. If one defines 
Zd as the value Mx/Tjj such that the X life time coincides with the age of the Universe 
(r = tu(zo)) then, for K -C 1, one has Zd — yJljK. Thus for K > 1 the X's will decay for 
Td = Mx/ Zd > Mx/\/2 and the inverse decays have to be take n into account . The kinetic 
equations (JSJ) and © are then generalized in the following wayl ' ' ' ' ' 



^ = -D(z)N x (z)+D(z)N^ 
dz 



D(z)Nx(z)+D(z)N^(z) (16) 



dN B-L ~ dN x , , , . , . 

~^z~ = -z—- W ^) N b-l{z)- (17) 

In the equation for Nx the second term accounts for the inverse decays that, remarkably, 
can now produce the X's. On the other hand one can see that a new term appears in 
the second equation for the asymmetry too, a wash-out term that tends to destroy what is 
generated from the decays. This term is controlled by the (inverse decays) wash-out factor 
given by 

m N eq 
W ID = -D-^^K, (18) 

z JV fe,Z 



'The equations <I16H and 1171 1 are actually not only accounting for decays and inverse decays but also for 
the real intermediate state contribution from 2 <-» 2 scattering processes. This term exactly cancels a CP 
non conserving term from inverse decays that would otherwise lead to an un-physical asymmetry generation 
in thermal equilibrium . 



where m is the number of baryons or leptons in the X decay final state (m = 1 in the case of 
leptogenesis) . Note that the decay parameter K is still the only parameter in the equations 
and thus the solutions will still depend only on K. They can be again worked out in an 
integral form . In the case of the B — L asymmetry one can write the final asymmetry as 

Ati Ari ~ L°° dz'W ID (z') ~ , s 

N B - L = N B ^ L e Jz > +eK{, (19) 

where now the efficiency factor is given by the integral 



Kf(K, z\) = — / dz 



dNx 
dz' 



■J? dz"W ID (z") 



(20) 



In the limit K — > the out-of-equilibrium case is recovered. In general one can see that the 
wash-out has the positive effect to damp a pre-existing asymmetry but also the negative 
one to damp the same asymmetry generated from decays, thus reducing the efficiency of 
the mechanism. A quantitative analysis is crucial and it is very useful to discuss separately 
the regime of strong wash out for K > 1 and the regime of weak wash-out for K < 1. 

2.3 Strong wash-out regime 

The strong wash-out regime is characterized by the existence, for K > 3, of an interval 
[z; n ,2; ou t] such that Wid ^ 1 and thus such that inverse decays are in equilibrium. Prac- 
tically all the asymmetry produced at z > z ou t is washed-out including, remarkably, an 
initial one. Moreover the calculation of the residual asymmetry is made very simple by the 
possibility to use the close equilibrium approximation given by 

dN x dN e * 2 Tr , , 

-sr ! *-iJi- = -ir, Wn >to- (21) 

In this way the integral in the Eq. I)20|) can be easily evaluated^ ' I Indeed, this can be 
regarded as a Laplace integral, that means an integral of the form 

dz'e-^ z '' z \ (22) 

o 

that receives a dominant contribution only from a small interval centered around a special 

value zb such that dip/dz = 0. In this way one can use the approximation of replacing 

Wid(z") with Wid(z") zb/z" in the Eq. 1)20(1 . With this approximation and assuming 

ZiCl the integral can be easily solved obtaining 

kfx—^fl-e—T*). (23) 

m K zb V / 



For large K 3> 1 and for m = 2 this expression coincides with that one found in 
The calculation of the important quantity zb proceeds from its definition, (dip/dz) ZB = 0, 
approximately equivalent to the equation 

W ID (z B ) = (-) (z B ) - - • (24) 
\ 7 / zb 



c Note however that the definition for K has to be used instead of K = (1/2) (r_D/-ff) z =i. 
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Figure 2: final efficiency factor as a function of the decay parameter K for thermal and dynamical initial X 

abundance. The solid lines are the numerical solutions of the Eq.'s (llfcil and 1171 . the short-dashed lines are 

the analytic results (cf. 12311 and I321 + 13311 ), the dot-dashed line is the power law fit Eq. 1541 . The dashed 

box is the range of values for K favored by neutrino mixing data in the case of leptogenesis (cf. 153H '). 



This is a transcendental algebraic equation and thus one cannot find an exact analytic 
solution (see'**' for an approximate procedure). However the expression 



z B (K) ~ 1 + 4.5 (miTr id e 



(25) 



provides quite a good fit that can be plugged into the Eq. (|2*5|) thus getting an analytic 
expression for the efficiency factor. At vary large K this behaves as a power law Kf oc i<f ~ 113 . 
In Fig. 2 we compare the analytic solution for K{ (cf. I|23|l ) with the numerical solution (for 
m = 1). One can see how for K > 4 the agreement is quite good. Note that the Eq. (|24|) 
implies that for large values of K one has zb — z on t, that particular value of z corresponding 
to the last moment when inverse decays are in equilibrium (Win > 1). In this way almost 
all the asymmetry produced for z < zb is washed-out and most of the surviving asymmetry 
is produced in the period just around the inverse decays freeze out, simply because the X 
abundance gets rapidly Boltzmann suppressed. An example of this picture is illustrated in 
Fig. 3 for K = 100 (from'**!). Instead of the abundance we plotted the deviation from the 
equilibrium value, the quantity A = Nx — N^. The deviation grows until the X's decay 
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Figure 3: comparison between analytical (short-dashed lines) and numerical (solid lines) results in 
the case of strong wash-out (K = 100) for |ei| = 0.75 x 10~ 6 . 



at z ~ Zd, when it reaches a maximum, and decrease afterwards when the abundance stays 
close to thermal equilibrium. Correspondingly the asymmetry grows for z < Zd, reaching 
a maximum around z ~ 1, and then it is washed-out until it freezes at zb — z ont . The 
evolution of the asymmetry Nb-l(z) can induce the wrong impression that the residual 
asymmetry is some fraction of what was generated at z — 1 and that one cannot relax 
the assumption z, < 1 without reducing considerably the final value of the asymmetry. 
Actually what is produced is also very quickly destroyed. A plot of the quantity ip(z, oo), 
as defined in the Eq. (|2'2j) and shown in Fig 4 (from^, enlightens some interesting aspects. 
This is the final asymmetry that was produced in a infinitesimal interval around z. It is 
evident how just the asymmetry that was produced around zb survives and, for this reason, 
the temperature Tg = Mx/zb can be rightly identified as the temperature of baryogenesis 
for these models. It also means that in the strong wash out regime the final asymmetry was 
produced when the X particles were fully non relativistic implying that the simple kinetic 
equations ()16|) and (|17|). employing the Boltzmann approximation, give actually accurate 
results and corrections from use of the exact quantum statistics can be safely neglected. 

This is not the only nice feature of the strong wash-out regime. Since any asymmetry 
generated for z < zb gets washed-out, one can also rightly neglect any pre-existing initial 
asymmetry N B _ L . At the same time the final asymmetry does not depend on the initial 
X abundance. In Fig. 5 we show how even starting from a zero abundance, the X's are 
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Figure 4: the function ip(z', 00 ) for if = 10 (solid line) and if = 100 (dashed line). 

rapidly produced by inverse decays in a way tha t w ell before zb the number of decaying 
neutrinos is always equal to the thermal number^". The final asymmetry does not even 
depend on the initial temperature as far as this is higher than ~ Tg and thus if one relaxes 
the assumption z% <C 1 to Z{ < zb — Azb, the final efficiency factor gets just slightly reduced 
(for example for Azb — 2 this is reduced approximately by 10%). 

Summarizing we can say that that in the strong wash out regime the reduced efficiency is 
compensated by the remarkable fact that, for 7] > Tg, the final asymmetry does not depend 
on the initial conditions and all non relativistic approximations work very well. These 
conclusions change quite drastically in the weak wash-out regime. 



2.4 Weak wash-out regime 

For K < 1 one can see that zb rapidly tends to unity (cf. (1251)). In Fig. 2 the analytic 
solution for the efficiency factor, Eq. Q23|) . is compared with the numerical solution. It 
can appear surprising that, in the case of an initial thermal abundance, the agreement 
is excellent not only at large K > 4, but also at small K < 0.4, with some appreciable 
deviation just around 0.4 < K < 4. The reason is that when the wash-out processes get 
frozen, the efficiency factor depends only on the initial number of neutrinos and not on 
its derivative and thus the approximation Eq. I|21|) introduces a sensible error only in the 
transition regime K ~ 1. 

The Eq. (|23j) can be easily generalized to any value of the initial abundance until one 
can neglect the X's produced by inverse decays. More generally, one has to calculate such 
a contribution and it is convenient to consider the limit case of a zero initial X abundance. 




Figure 5: fast thcrmalization of the X abundance in the strong wash-out regime. The final Nb-l 

abundance (for E\ — 0.75 x 10 -6 ) is the same as in the case of an initial thermal abundance (cf. Fig. 

3) and it is independent on the evolution at z -C Zb- 

The X production lasts until z = z cq , when the abundance is equal to the equilibrium value, 
such that 

N x {z cq ) = N x i (z eq ). (26) 

At this time the number of decays equals the number of inverse decays. For z < z eq decays 
can be neglected and the Eq. (|TB|) becomes 

dN x 



dz 



D(z)N e x \z). 



For z«l one then simply finds 



Nx(z) 



K 

~6 



z 3 . 



(27) 



(28) 



In the weak wash out regime the equilibrium is reached very late, when neutrinos are already 
non relativistic and z cq S> 1. In this way one can see that the number of Nx reaches, at 
z ~ z cq , a maximum value given by 



N x ( Zeq )-N(K) = ^K. 



(29) 



It is possible to interpolate between the two asymptotical regimes getting a global solution 
for any z < z eq . For z > z cq inverse decays can be neglected and the X's decay out of 
equilibrium in a way that 



N x (z > z eq ) ~ N x (z cq ) e 



- r dz'D(z') 



(30) 



Let us now consider the evolution of the asymmetry calculating the efficiency factor. Its 
value can be conveniently decomposed as the sum of two contributions, a negative one, k7 , 
generated at z < z eq , and a positive one, nf , generated at z > z eq . In the limit of zero 
wash-out we know that the final efficiency factor must vanish, since we are assuming an 
initial zero abundance. This implies that the negative and the positive contributions have 
to cancel each other. The effect of wash-out is to suppress the negative contribution more 
than the positive one, in a way that the cancellation is only partial. In the weak wash-out 
regime it is possible in first approximation to neglect completely the wash-out at z > z ccl . 
In this way it is easy to derive from the Eq. (j20|) the following expression for the final 
efficiency factor: 

Kt =* N{K) - 2 ( 1 - e ~h N ( K )) . (31) 



One can see how it vanishes at the first order in N(K) oc K and only at the second order 
one gets k f ~ (9vr 2 /64) jf 2 H3 

2. 5 Final efficiency factor: summary 

Generalizing the procedure seen for the strong wash-out it is possible to find a global solution 
for Kf(K) valid for any K. The calculation proceeds separately for k~ and k + and the final 
results are given by 

k7(K) = -2 e-^W ( e pT(*) _ i) ( 32 ) 

and 

The function N{K) extends, approximately, the definition of N(K) to any value of K 

N(K) =: ; N( ^L 2 ■ (34) 

The sum of the Eq.'s (|33|) and (|32[) is plotted, for m = 1, in Fig. 2 (short-dashed line) and 
compared with the numerical solution (solid line). 

We can now outline some conclusions about a comparison between the weak and the 
strong wash-out regimes. A large efficiency in the weak wash-out regime relies on some 
unspecified mechanism that should have produced a large (thermal or non thermal) X 
abundance before their decays. On the other hand the decrease of the efficiency at large 
K in the strong wash-out regime is only (approximately) linear and not exponential!^. 
This means that for moderately large values of if a small loss in the efficiency would be 
compensated by a full thermal description such that the predicted asymmetry does not 
depend on the initial conditions, a nice situation that resembles closely the Standard Big 
Bang Nucleosynthesis scenario for the calculation of the primordial nuclear abundances. 



3 Leptogenesis 

Let us see now how the results that hold for generi cal baryogenesis models from heavy particle 
decays get specialized in the case of leptogenesis'^. This is the cosmological consequence 
of the seesaw mechanism, explaining the lightness of the ordinary neutrinos through the 
existence of three new heavy RH neutrinos N\,N2, N3 with masses respectively Mi < M2 < 
M3, much larger than the electroweak scale. The simple seesaw formula, 

1 T 

m v = - m Djrm D , (35) 

relates the neutrino mixing matrix m v to the RH neutrino mass matrix M and to the Dirac 
neutrino mass matrix mp = hv generated by the Yukawa coupling matrix h, where v is 
the Higgs vacuum expectation value. Both light and heavy neutrinos are predicted to be 
Majorana neutrinos. All mass matrices are in general complex and this provides a natural 
source for the CP asymmetry while the new RH neutrinos are the natural candidates to 
play the role of the X particles. In this case things are apparently more complicated since 
there are three of them. We will assume that the decays and inverse decays of the two 
heavier neutrino decays do not influence the value of the final asymmetry. This assumption 
holds for example either if the asymmetry produced by the two heavier RH neutrinos is 
negligible or if this is produced and then washed out by the inverse decays of the lightest 
(heavy RH) . In this way we can straightforwardly apply the general picture of baryogenesis 
from X decays to leptogenesis, with the N\s playing the role of the X particles. 

3.1 Decay parameter and neutrino masses 
The total N± decay width is given by 

t _ miM[ 

where the effective neutrino mass is defined as^^ 

(TrJ D m D )n ,„„, 

mi = Ml ■ (37) 

It is then easy to see that the decay parameter is related to mi'™ by the following relation 

K = ^, (38) 

where the equilibrium neutrino mass m* can be written as 

v 2 
m ± = ~ 10" 3 eV , (39) 



with the quantity M* given by 



3^ Mpi _ ., , (i| 



M * = 77^72 ^=- 3xl0GeV - ( 4 °) 



It is quite non trivial that the value of m* is close to the neutrino mixing mass scales and 
we will show soon the relevance of this result. For the moment note that the value of 
m* is independent on the well known success of the seesaw mechanism in explaining the 
atmospheric and solar neutrino mass scales and this is why we wrote m* in a sort of seesaw- 
like form, introducing the scale M*. Apart from the very generical consideration that the 
logarithm of M* is expected to be close to the Planck scale, this is not related to the grand 
unified scale, rather to the expansion rate at the baryogenesis time d . 

Let us now assume that the simple decays plus inverse decays picture studied in the 
previous section is a good approximation of leptogenesis. It is then crucial to determine the 
value of the the effective neutrino mass fh\, and thus, from the Eq. (|38[) . the value of the 
decay parameter K, in order to answer the important question whether leptogenesis lies in 
the strong or in the weak wash-out regime. 

It is always possible to work in a basis in which the heavy neutrino mass matrix is 
diagonal, such that M = diag(Mi, M2, M3) = Dm- Moreover one can also simultaneously 
diagonalize the light neutrino mass matrix m u by mean of the unitary MNS matrix U, such 
that 

[/t mv U* = -D m . (41) 

In this way the seesaw formula Q35JI gets specialized in the following way: 

D m = U^m D D M l m T D U*. (42) 

This expression can be also re-casted as an orthogonality condition, 

nn T = tt T Q = I, (43) 

for the $7 matrix defined as'==' 

Sl = D m l l 2 tfm D D- M l/2 (44) 

and whose matrix elements are then simply given by 

n« = 4% ' (45) 

where h = U* h. The Q matrix is fully determined by three complex parameters. Four 
of them are needed to fix the three first column entries f^n, f^l an d ^31, particularly 
important for leptogenesis. This because if one inverts the relation (|44jl. in a way to get an 
expression of mo in term s of £1, and plugs it into the effective neutrino mass definition (cf. 
(|37|)h then one easily gets'**" 

?h\ = Till \£ln\ + 7TJ2 |^2ll + rn 3 l^3ll ■ (46) 



It is then quite curious that the value of M* is just the value of the supersymmetric unification scale. 



From the orthogonality of 0, it follows that mi > m\. This is the only fully model indepen- 
dent restriction on fh\. For configurations such that 



Ei" 2 



'in 



£<* 



(47) 



one has rh\ < 711,3. Models with rh\ 3> 7713 rely on the possibility of strong phase cancella- 
tions. 

Neutrino mixing data provide two important pieces of information on the neutrino 
mass spectrum. In the case of normal hierarchy one has m| — m 2 . = Am atm and m 2 , ~ fn\ = 
Am 2 ol . In the case of inverted hierarchy m| — m 2 . = Am 2 ol and m 2 , — m 2 = Am atm . The 
third, still undetermined, independent information, the absolute neutrino mass scale, can be 
conveniently expressed in terms of the lightest neutrino mass mi . The two heavier neutrino 
masses are then given, for normal (inverted) hierarchy, by 

ml = m\ + ml %m , (48) 

m 2 = m 2 + Am 2 ol (Am 2 tm ), (49) 



where we defined m atm = JAm^ + Ara s 2 ol . The latest measurements give^^ 

Aml tm = (2.6 ± 0.4) x 10" 3 eV 2 , (50) 



usi 



and for solar neutrinos 



from which it follows that 



Am s ^~(7.11^xlO- 5 )eV^, (51) 



m atm = (0.051 ± 0.004) eV . (52) 

These relations imply that for m 2 3> m 2 tm neutrinos are quasi-degenerate (171,3 — m 2 — Trii), 
whereas for m 2 <C m atm they are hierarchical (ttii <C 7112,777,3). 

For fully hierarchical neutrinos (mi = 0) there is practically no restriction on rh\. 
However the case rh\ <C 777,2,771,3 requires l^il << 1 anc ^ l^iil << vtiijmz- This situation 
cannot be excluded!*^ but, because of the observed large mixing angles in the mixing matrix 
U, it relies on a fine tuning between the U and mo matrix elements (cf. (|44|)). such that 
the off-diagonal terms are very small. This qualitative and general argument is supported 
by different investigations on specific models or classes of models for which typically one 
finds m so \ ~ ?ii2 ^ fh\ < 7713 ~ Tn a t m '=^. Therefore, in the case of normal hierarchy one 
has that the favored range for the mi value is given by C(m so i) < fh\ < 0(m a t m )> that in 
terms of the decay parameter (cf. (J3H|)) gets translated into the range 

0(K sol ~ 7) < if < 0(K atm ~ 50) , (53) 

while for inverted hierarchy the situation is even simpler since rh\ = 0(m atm ) and K = 
0(K a t m ~ 50). One thus arrives to the interesting conclusion that neutrino mixing data 



favor leptogenesis to lie in a mildly strong wash out regime, strong enough to benefit from 
the advantages we discussed, independence on the initial conditions plus minimal theoretical 
uncertainties, but not too much to result in an untenable efficiency loss. This conclusion 
derives because both the two independent experimental quantities, m so \ and m a t m > are about 
ten times m+ and so now one can better appreciate the nice matching of the theoretical 
quantity m* with the experimental data 6 . In the range K so \ < K < if atm a good fit of the 
final efficiency factor (cf. Eq. I|23|)) is given by the power law 

shown in Fig. 2 (dot-dashed line). These conclusions hold under the assumption that 
leptogenesis is well approximated by the simple decays plus inverse decays picture and we 
have now to verify whether they are drastically modified or just corrected by the account 
of N\ scatterings and AL = 2 processes. 

3. 2 Scatterings 

The iVi's can also be destroyed or produced in AL = 1 scatterings involving the top quark. 
These are mediated by the Higgs and can occur in the s channel, like N\ + I <-> t + q, or 
in the t channel, like Ni + t <-> I + q. The account of these processes modify the kinetic 
equations lfT5|l and (|T7|) in the following way: 

^ = -(D + S)(N Nl -N%), (55) 

^± = - £l D {N Nl - N%) - (W ID + W AL+1 ) N B ^ L . (56) 

Note that scatterings have two effects: they contribute both to the neutrino production 
(the S function) and to the wash-out (the Wal+i function). 

The first one is important in the weak wash-out regime. As one can see from the 
Eq.'s (|55|) and (|56|) . the production of the JVi's from the S function is not associated to a 
production of the asymmetry, simply because these processes do not violate CP. In Fig. 5 
(froml"!) we show an example of N\ production for fh\ = 10 -5 eV (K ~ 0.01), comparing 
the case when scatterings are included with the case when they are neglected. It can be 
seen how at z = z eq the number of neutrinos is approximately doubled while the final 
asymmetry is two orders of magnitude larger. The reason is that the neutrino production 
from the scatterings is not associated to a production of a negative asymmetry. On the other 
hand all produced neutrinos yield a positive contribution when they decay. The expression 
(|Hljl for the final efficiency factor in the weak wash-out regime gets thus modified in the 
following way at the first order in K 

N Nl D 



Kf 



D + S 



-N(K)ocK. (57) 

e Note that this is also a consequence of the recent exclusion of the low solution in the solar neutrino data, 
that would have implied K so \ <^. 1. 




z=M/T 

Figure 6: comparison between the case when scatterings are included (thick lines) with the decays 
plus inverse decays picture (solid thin lines) for rh\ = 10 -5 eV. 



If scatterings are switched off the negative and the positive contribution cancel at the first 
order like we saw already. If S 7^ the positive term is enhanced while the negative one 
remains unchanged and in this way the sum does not vanish any more. Hence this effect 
makes more efficient the asymmetry production at small K, without having to assume an 
initial thermal abundance. There is however a drawback. The final result is quite sensitive 
to the theoretical assumptions. The scattering cross section depends on the ratio, M^/Mi, of 
the Higgs mass to the RH neutrino mass. The case depicted in Fig. 5 is for Af^/Mi = 10 -5 . 
For smaller values of th is ratio the result does not change much. However it has been 
recently pointed out!" that the Higgs mass is better described by its thermal mass such 
that Mh/Mi ~ 0.4/2. The relevant values of z for neutrino production are z — 1 and so the 
ratio Mfr/Mi ~ 0.4. Such an high value has the effect to suppress heavily the S term and the 
suppression is made even stronger by the account of the running of the top Yukawa coupling 
at high temperature. In this way the simple decays plus inverse decays picture is practically 
recovered. On the other hand inl " I ' 1ft has been noticed how scatterings involving gauge 
bosons should also be included. These scatterings yield an additional contribution to the S 
function such that the final result is between a situation where scatterings are neglected and 
one where scatterings involving top quark and small M^/Mi are taken into account. The 
conclusion is that in the weak-wash out regime the theoretical uncertainties are such that 
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Figure 7: efficiency factor when scatterings are included. 
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it seems that any result between the simple decays plus inverse decays picture, for which 
Kf oc K 2 or a behavior K{ oc K (cf. Eq. (|57|) ) cannot be firmly excluded at the moment. 
These large theoretical uncertainties, represented in Fig. 6 with the short-dashed region, 
are in addition to the model dependence in the description of the initial conditions. 

In the strong wash-out regime all difficulties get considerably reduced. The theoretical 
uncertainties in the description od scatterings can change the final efficiency factor no more 
than (20 -r 30)% and this is clearly shown in Fig. 6, where at large K the (thin solid line) 
range shrinks considerably compared to the (short-dashed line) range at small K. This 
because the thermal abundance limit is saturated at z eq <C 1 anyway and therefore the 
number of decaying neutrinos does not depend on the S function. A residual source of 
uncertainty is still present because of the scattering contribution, Wal=i, to the wash-out. 
The effect of this term is however small, for the simple reason that in the strong wash 
out regime the surviving asymmetry is produced sharply around z^ > 1 and at such low 
temperatures inverse decays are dominant compared to scatterings. The conclusion is that 
in the strong wash-out regime the simple decays plus inverse decays picture does not get 
modified by scatterings within the theoretical uncertainties. It has been also pointed out'***' 
that an accurate description of the dynamics of sphalerons in converting the lepton number 
into a baryon number is expected to lead to a suppression of the final asymmetry of a 
0(1) factor and since this is currently neglected it gives an additional contribution to the 
theoretical uncertainties. Taking into account all these effects, an expression for the final 
efficiency factor in the strong wash out regime that accounts for the theoretical uncertainties 



is given by the power law' 



.3 



9 /l(T 2 eV\ 
Kf = (2 ± 1) 10~ 2 — = . (58) 



The central value corresponds to the curve represented in Fig. 5 with circles (more precisely 
this is obtained for a power law mj" ' ), while the range that is spanned by the error, 
corresponds approximately to the thin solid line area. The upper values of this range is the 
power law (J5I|) , that well describes the simple decays plus inverse decays picture where the 
wash-out from scatterings is neglected-'. 

3.3 AL = 2 processes 

There is another important contribution to the wash-out term arising from the AL = 2 
processes like IH <-> IH and med iate d by the RH neutrinos. In the non relativistic regime 
this contribution tends simply to'""' 

uj ( M\ \ ( ' m \ 



with m 2 = mf+m^+m^, and dominates on the other Boltzmann suppressed wash-out terms 
arising from inverse decays and scatterings. A well known problem is that at temperatures 
T ~ Mi one has to be sure that the cross section of AL = 2 processes does not double 
count the on-shell contribution already accounted by inverse decays followed by decays to a 
final state with opposite lepton number (i.e. I + H —> Ni —> 1 + H). In ' it has been found 
that the subtraction procedure usually employed in the previous literature gives arise to a 
washout AL = 2 term that is very well approximated by the asymptotical non-relativistic 
limit Eq. (|59|) plus a term that is just half the washout from inverse decays. In^J this 
second term has been shown to be spurious and to disappear when a proper subtraction 
procedure is employed. This result has been confirmed in'"'. Therefore the effect of the 
AL = 2 processes is entirely well approximated by its non-relativistic limit. It is easy to 
see that for Mi <C 10 14 GeV (0.05 eV/m) 2 , this term can be neglected. Thus, for sufficiently 
small neutrino masses and in the strong wash-out regime, we can conclude that leptogenesis 
is well approximated by a simple decays plus inverse decays picture. 

3.4 CP asymmetry and seesaw geometry 

So far we concentrated on the kinetic theory of leptogenesis and we have seen how neutrino 
mixing data favor a very simple regime in which predictions are model independent and the- 
oretical uncertainties are minimized. We have now to answer the crucial question whether 
the resulting final asymmetry can explain the measured CMB value (cf. (|Tjl). The thermo- 
dynamical point of view, i.e. the efficiency factor, is not enough to answer this question, 

■''The result obtained in'" corresponds to this situation because at 2 ~ 2b > 1 the Higgs thermal mass 
suppresses the wash-out from scatterings involving the top quark while the contribution from scatterings 
involving gauge bosons is negligible. 



since one needs to know the value of the CP asymmetry too. This is a specific leptogenesis 
issue that concerns what can be called the seesaw geometry. 

A perturbative calculation from the interference between tree level and vertex plus self 
energy one-loop diagrams yields'^ 
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The function fy, describing the vertex contribution, is given by 

1 + x 



fv(x) 



1 - (1 +x) In 



while the function fg, describing the self-energy contribution, is given by 

fs(x) = 



x 



(60) 



(61) 



(62) 



In the limit i>1, corresponding to have a mild RH neutrinos mass hierarchy with M 2 ^ 3> 
Ml, one has 

Mi)+/s(i)=!~. (63) 

In this limit and barring strong phase cancellations'"^ the expression (|6U|) simplifies into'™ 

3 Im [(rf h M- 1 h T h*) u ] 



£l ^~ 16vr (/i/it)n 

Replacing h with fl (cf. l)I5|) ) one then gets^^ 

3 Mi m atm 
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where we introduced the convenient quantity 

(3(m 1 ,m 1 ,n 2 jl ) 



-Tr. o f3(m 1 ,m 1 ,n jl ) 



iriatm J2j m j l^jll 



(64) 



(65) 



(66) 



The final asymmetry is proportional to the product of the CP asymmetry times the final 
efficiency factor that, in the simplified decays plus inverse decays picture, depends only 
on the effective neutrino mass fh\. The expression (|65|) shows that the CP asymmetry 
depends on the three complex numbers f^ and thus it introduces a model dependence 
in the prediction of the final asymmetry that one was hoping to have removed in the 
calculation of the final efficiency factor. It is however possible to maximize the absolute 
value of the CP asymmetry respect to the 'geometrical' parameters il^ 1; thus finding a 
non trivial maximum ^^(Mi, mi, mi) depending only on M\, rh\ and mi. One can then 
define an effective leptogenesis phase 5l such that the expression ()66|) can be re-casted in 
the following way 



(3(m 1 ,mi,VL 2 jl ) = (3v 



.(mi, mi) sm5 L (m 1 ,mi,£l'j 1 ) 



(67) 
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Figure 8: Seesaw geometry. Different configurations of 0,% (see text for explanation) with the same 

value of rh\ . 



The maximum of the absolute value of the CP asymmetry and of the function (3 are thus 
realized for those particular geometrical configurations, corresponding to some O^'s values, 
s uch that sin<5k = 1. A general procedure for the calculation of e™ ax and svo.5l is presented 
i n l"i2l Here we just sketch some general features and describe two particularly interesting 
limit cases. 

If one represents the three f& in the complex plane, the orthogonality condition fixes 
the sum of the three to start from the origin and to end up onto the real axis at the point 
ReQ^ a ^?i) = 1) as shown in Fig. 8 for a generic configuration (solid line arrows). Using 

Xj + iYj and using the definition of rh\ (cf. 



the orthogonality condition, defining Q, 2 ^ 
EHJ, this can be re-casted as 



Am| 2 Y 3 + Am? i Y 



m at m mi 

with Y = Y2 + 13. The absolute value of (3 has to be maximized for rh\ constant. In general 
one always finds that """^ 

Anax(mi,mi) = /5 max (mi) /(mi, mi) < 1, (69) 

with /3 max (mi) = (m 3 - mi)/m atm , /(mi, mi) < 1 and /(mi,oo) = 1. 

An interesting limit case is that of /u% hierarchical neutrinos for mi = 0. In this case 
Anax = 1 an d there is no global suppression. Moreover one has mi = mi |f^lil + 7TJ -3 l^iil an d 
thus, for any change configuration such that |f2|il and l^f^l are constant, the quantity rh\ is 
also constant while If^l can be arbitrarily modified. Hence |/3(mi,mi,f&)| is maximized 



for configurations such that X2 = X3 = 0. It is then easy to see that it is further maximized 
for Y2 = and Y3 = fhi/m^, corresponding to the configuration shown in Fig. 8 with dotted 
line arrows. In this case one has very simply /(0, mi) = 1. Therefore, the case m i = 
corresponds, for a fixed M±, to an absolute maximum of the CP asymmetry given byl ' ' 
(cf. dHSJ) and HU)) 

ri-W) = -A- Mimatm ~ 10- 6 f Ml "\ ( matm 1 (70) 

1 l lj 16vr w 2 V 10l ° GeV y V0.05eV; { } 

Note that with this last definition of e™ ax (Mi), together with the expressions 1)6 7|) and 1)69)1 . 
the Eq. 1)65)1 for the CP asymmetry can be re-casted like 

£l = -ei MX (A''i)Anax(mi)/(mi,mi) sin^fm!,^,^) , (71) 

showing the sequence of different maximization steps. 

In the quasi- degenerate limit the expression (|46)1 for m-i becomes simply m-i ~ mi ^ ■ |f2?i| 
Thus the condition mi = const is equivalent to select all those configurations for which 
^\- \Q 2 ji\ is constant. Hence it is straightforward to conclude that \/3\ is maximum for a 
configuration such that Y2 = -X2 = and X3 = 1/2, shown in Fig. 8 with dashed line 
arrows. Using that in the quasi degenerate limit m 2 tm ~ 2 mi (m.3 — mi), one obtains 9 



/(mi,mi) = Jl-^i. (72) 

V m i 

Note that in both the two limit cases the maximum CP asymmetry is obtained for config- 
urations su ch t hat O21 = 0. It is possible to show that this result holds in general, for any 
value of mi™ 

Therefore, for maximal CP asymmetry (sinJ^ = 1), one can still express all predictions 
in terms just of M\,fh\ and mi. In particular it is possible to express the CMB constraints, 
for all neutrino models, just in terms of these three parameters. For specific models it can 
happen of course that s\n.5i < 1 and the constraints get, in general, more restrictive. 

3.5 CMB bound 

From the Eq.'s (jTHj) and 1)19)1 one obtains for the predicted baryon to photon number ratio 

r/B =de\K{, (73) 

where the quantity d is defined as 

d=^. (74) 

7 

In the Standard Model case one has a sp h = 28/79, while the number of photons at recom- 
bination, assuming a standard thermal history, is given by 

™ = i^SM = 4697 36 

Z_ 3 5rec 129 V ; 

9 This limit expression has been first shown in"" using the approximation m so i = 0. Here we derived it 
in a more general way. 



and thus d~ 0.97 x 1CT 2 . 

The maximum baryon asymmetry T] , ^ ax (Mi,fhi,mi) is defined like the asymmetry cor- 
responding to the maximum CP asymmetry 

ef^Mi.mi.mi) =eH«i)^axK)/(m 1) m 1 ). (76) 

The CMB bound is then simply equivalent to require 

CfMbmLm^^ (77) 

and therefore will yield constraints on the space of the three parameters M\, m\ and rh\. 

3.6 Lower bounds on the lightest RH neutrino mass and on the reheating temperature 
We have seen that the absolute maximum of the CP asymmetry is obtained for mi = 0. For 

or 

mi > the function /3 max (mi) suppresses the CP asymmetry'^'. Furthermore the AL = 2 
wash-out term gets enhanced when the absolute neutrino mass scale increases (cf. (|59|0 , 
Therefore, the maximum baryon asymmetry r/^ ax is maximal when mi = 0. In this case the 
allowed region in the space of the parameters M\ and fh\ and compatible with the CMB 
constraint is maximum^ and one finds an interesting lower bound on the M\ value' ' ' 
just plugging the expression ((70)) into the CMB constraint (cf. (J77J)) 

Ml > M r - 116 ' "' ri%MB 



d 3 m atm «f 

7?g MB \ / 0.05eV 
6xl0- 10 ) I m atm 



- 6.4xlO«GeV(^^)(^^)^. (78) 



For an initial thermal abundance and in the limit m-i/m* — ► 0, one has, by definition, Kf = 1 
and so one finds 

Mi > (6.6 ± 0.8) x 10 8 GeV > 4 x 10 8 GeV , (79) 

where the last inequality is the 3<r bound and we have used the experimental values Eq. 
(PQ) and Eq. (|52|). The case of a dynamically generated N\ abundance is more significative 
and in this case the peak value Kf ~ 0.18, implies 

Mi > (3.6 ± 0.4) x 10 9 GeV > 2 x 10 9 GeV . (80) 

The most interesting situation corresponds to the range (m so \ H-m atm ), for which the power 
law Eq. (|58~|) can be used for Kf, thus giving 

/ ~ \ 1.1(1.2) 

Mi > [3.3(2.2)±0.4(0.3)]xl0 10 GeV f ^z^J > [1.5(1)-10(9)] x 10 10 GeV , (81) 

where we have used the central value (upper value) in the Eq. (|58[), The Mi lower bound 
can be translated into a lower bound on the initial temperature T\ that, within inflationary 
models, can be identified with the reheating temperature, corresponding to that temperature 



below which a radiation dominated regime holds'^. So far we assumed that this is much 
larger than M±. If one relaxes this assumption then the final efficiency factor gets reduced. 
For small rh\ the threshold value is given approximately by M\ itself since below this 
temperature either, assuming an initial thermal abundance, the Aq abundance is thermally 
suppressed or the production gets considerably suppressed for an initial zero abundance. 
Therefore, for small values rh\ < 10 -3 eV, the same bounds (f?9"|) and (|%U|) apply also to the 
reheating temperature. 

In the more interesting case of strong wash-out, since the 90% of the surviving abundance 
is produced in an interval z ~ zb ± 2, then the reheating temperature can be ~ zb — 2 
times lower than Mi™ without any appreciable change in the final predicted asymmetry. 
In the interesting range (m so i,m a tm) one has that zb spans between 6 and 8 and thus the 
bound Eq. (|81|) gets relaxed from 4 to 6 times giving 

Ti > [4 (2.5) x 10 9 - 2 (1.5) x 10 10 ] GeV . (82) 

This is another interesting result showing how the support of neutrino mixing data to the 
range rh\ ~ 0(m ao i, m a tm) not only makes leptogenesis working in a simple and predictive 
way but also how the loss in the efficiency is compensated by a non relativistic production 
of the final asymmetry such that the lower bound on the reheating temperature gets just 
slightly more restrictive compared to the small fh\ range. Note that in case of modifications 
of the theoretical assumptions such that the maximum baryon asymmetry rj r ^ a,x — > £,r/ B ia,x , 
one has correspondingly Mf m ,T™' m -> M™ in /(,,T t min /£. 

3.7 Upper bound on the absolute neutrino mass scale 

For large values of the absolute neutrino mass scale the AW wash-out term cannot be 
neglected. The final efficiency factor can be calculated in the approximation that AW starts 
to be effective for z > zb, when the asymmetry generation from decays already stopped. 
This is a very good approximation in the strong wash-out regime and since rh\ > m\ it 
does not introduce any restriction for m\ > m* ~ 10~ 3 eV. Within this approximation one 
has simply 

«f(mi, Mirrr) = «f(mi)e z b \w^~Giv) Vcv; ; ( 83 ) 

where Kf(fhi) is the efficiency factor calculated in the regime of small neutrino masses 
neglecting the AW term. We use the simple limit kf = 2/(zb K) that corresponds to 
neglect, conservatively, the contribution of scatterings to the wash-out and approximately 
to the upper values in the Eq. I|58|). A search of the peak values (Ml, mi) for which the 
maximum baryon asymmetry r/g ax has an absolute maximum yields 



V CMB ~ 37/2 m 4 



(84) 



with the constant \ ~ 1.6 eV 3 . Hence t he CM B bound implies an interesting constraint on 
the neutrino mass, given by mi < 0.1 eV* ' I A precise calculation has to take into account 



the running of neutrino masses'^. The atmospheric neutrino mass scale at temperatures 
T ~ 10 13 GeV is higher than at zero temperature. On the other hand the bound on 
neutrino masses that is obtained at large temperatures gets lower when calculated at low 
temperatures. This second effect is dominant and thus the account of neutrino mass running 
will make the neutrino mass bound more restrictive. The smallest effect is obtained for an 
Higgs mass Mh — 150 GeV and makes the bound 20% more stringent. Taking into account 
this effect one then obtains the 3cr bound'"' 



m. 



<0.12eV£ 1/4 . (85) 



Note however that a (1 figure) bound rrii < 0.1 eV conservatively accounts for the theo- 
retical uncertainties h . As defined in the previous subsection, a value of £ ^ 1 describes a 
possible variation of the maximum baryon asymmetry in the case of modified theoretical 
assumptions, like for example in the supersymmetric case that will be studied in the next 
subsection, or in the presence of possible different effects, lik e an enhancement of the CP 
asymmet ry due to a degenerate heavy neutrino spectrum'^ that would relax the bounds 
' ^ l ' ^ ', or simply for a variation of the input values of the experimental quantities (note 
that £ ex fn'lx,ml' r l < B MB )- However the strong suppression of the baryon asymmetry for an 
increasing absolute neutrino mass scale (tjb oc 1/mf) makes the bound quite stable^. It is 
important to realize that the bound can be evaded but not trivially. This means that a mea- 
surement of a value of the absolute neutrino mass scale above the leptogenesis bound will 
necessarily imply some drastic modifications of the minimal leptogen esis sce nario. These 
include particular neutrino models within the simple seesaw formula' '', non thermal 
leptogenesis scenari os'^ or a non minimal seesaw formula, like that one arising in theories 
with a triplet Higgs . 

3.8 The supersymmetric case 

Leptogenesis can be also studied within the minimal supersymmetric standard model (MSSM) 
' ' ' '. In this case the asymmetry is generated not only from the N\ decays but also from 
the decays of their scalar partners, the iVf's and their antiparticles iVf , with the same mass 
M\. Since the decay width and thus also the inverse decay wash-out term are the same, 
these yield an additional equal contribution 1 . Therefore, from a thermodynamical point of 
view, the Aq's and the Aq's will play the same role and it is simply like if the 'X-abundance' 



h If instead of the upper values we were using the central (lower) values in the Eq. Q580 . then the bound 
would have been ~ 0.005 (0.01) eV more stringent. As we said already, these values could arise from a 
possible contribution to the wash-out from scatterings or for the account of spectator processes. Note also 
that we use d th e latest published results on Am atm ' — '• Preliminary results from the SuperKamiokande 
collaboration^] find AmL = (1.3 - 3.0) x 10" 3 eV 2 at 90% c.l. (best fit AmL = 2.0 x 10" 3 eV 2 ), from 
which, at la, m atm = (0.045 ± 0.006) eV, implying a ~ 0.01 eV more stringent bound rrii < 0.11 eV. IrJ— ■*& 
bound rrn < 0.15 eV, when this second value of m a tm is employed, has been obtained. A difference ~ 0.02 eV 
can be ascribed to a different estimation of the running of neutrino masses. The remaining 0.02 eV difference 
can be safely included in the account of the theoretical uncertainty. 

l The discussion in this subsection is fromM 



gets doubled. We will still track the B — L asymmetry in the co-moving volume containing, 
on average in ultra-relativistic thermal equilibrium, one RH neutrino N\ and thus now also 
one N\. In this way we will still write the final baryon asymmetry as in the Eq. (|13|) but 
now d = 2a sp } 1 /-/vT ec , with the additional factor 2 taking into account the contribution from 
the N^s 3 . Let us analyze how the different quantities involved in the calculation of the 
final asymmetry get modified from the SM to the MSSM case. To this aim it will prove 
convenient to introduce the variations 

XMSSM 

& - x SM ^ 

for any quantity X. The sphaleron conversion coefficient is given by o sp h = 8/23 and thus 
it is almost unchanged. The number of degrees of freedom is given by gMSSM = 915/4 and 
this approximately doubles the number of photons at recombination given by 
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In this way one gets d ~ 0.89 x 10 2 , almost unchanged compared to the SM case 

(£d ~ 0.92). The CP asymmetry in the MSSM case is given byl^ 



r2\ / -\/r2\-\ 



s^h (hWi x r [m!J +fs {m!J 



fv{x) = ^\n[ l -±^) (89) 



=2,3 

where 

and _ 

fs(x) = ^- . (90) 

1 — x 

In the limit i>1, corresponding to have a mild RH neutrinos mass hierarchy with Mf 3> 
Ml and barring strong phase cancellations, one has 

M*) + //?(*) — 4= ( 91 ) 

Jx 



and thus, compared to the SM case (cf. (|60|) and (|63|)). the absolute value of the CP 
asymmetry and hence also of its maximum e" ia,x get doubled (£ e = 2). 

Let us now calculate how the efficiency factor at small values of Mi m 2 gets modified 
in the MSSM case. We have seen that SM leptogenesis is well approximated, in the strong 
wash-out regime, by the simple decays plus inverse decays picture. It is then interesting 
to study how the results change within such a simple picture. Since the JVi's can now 
decay in two new channels (N\ — ► lh,tfh) that give exactly the same contribution to the 
decay width as the other two standard ones, this gets doubled compared to the standard 

-'The other possibility would have been to choose a halved co-moving volume in a way that the factor 2 
was absorbed in A^ ec . 



case (cf. (|36|l). This makes lifetime shorter and the inverse decays wash out rate stronger. 
However the increase of the degrees of freedom makes the expansion faster and this partially 
compensates. Recalling the definition of the equilibrium neutrino mass, Eq. (j39j) . one has 
simply 

^ D 

implying m MSSM ~ 0.8 x 10~ 3 eV. Therefore now one has that the transition to the strong 
wash-out regime occurs for slightly smaller values of rh\ (about l/\/2 smaller) and thus for 
the decay parameter one has £k — v2. The range of K values favored by neutrino mixing 
data will thus be given by 

C^MSSM „ 1Q) < K MSSM < (K^i SM ~ 65) , (93) 

to be compared with the Eq. Q53JI in the SM case. In the strong wash out regime the 
efficiency factor, calculated in the decays plus inverse decays picture, will be still given by 
the expression Eq. 1|23|) . For K in the range (|93|) one has K{ oc K~ 1A5 (cf. (|54|l) and thus 

^-e^-^=- (94) 

Assuming that the effect of scatterings in the MSSM goes in the same direction as in the 
SM and thus that the result ^ f . ~ 1/ ' \[2 holds approximately also when scatterings are 
included, one can write the analogous of the Eq. (j20[) in the MSSM case 

K MSSM = (L5 ± 0>7) 1Q -2 / l^V \ 

V mi y 

A more detailed analysis is needed to verify the role of scatterings in the MSSM. Note that 
again the upper values of this range correspond to the power law (|54(l translated in the 
MSSM case (i.e. by replacing m^ M -> m* ISSM ). 

Let us now investigate the consequences for the lower bound on Mi. The Eq. (|78|) will 
now become 

Mi > Mf 1111 = -!£ 

d 3 m atm K{ 

„ 3.5 x 10» GeV ( "§"",„) (^^) «{■' . (96) 

From this expression one obtains, for the case of initial thermal abundance, zero initial 
abundance and in the strong wash-out regime respectively, the constraints 

Mi > (3.7 ± 0.4) x 10 8 GeV > 2.5 x 10 8 GeV , (97) 

Mi > (1.7 ± 0.2) x 10 9 GeV > 1.1 x 10 9 GeV (98) 



and 

/ ~ \ 1.1(1-2) 
Mi > [2.4 (1.6) ± 0.2 (0.15)] x 10 10 GeV ( j^-^ J > [1.5 (0.9) - 10 (8)] x 10 10 GeV , 

(99) 
where we have used the central (upper) value in the Eq. 1)95)1 . These have to be compared 
with the constraints obtained in the SM case (cf. Eq.'s (|79jl. (|8()jl and (|81|l). In the first 
two cases the constraints are approximately twice looser, because of the CP asymmetry 
enhancement, while in the case of strong wash out one has £M min — (£e £«_*) — 0.7 for the 
central value, while the 3<t bound remains practically unchanged because the experimental 
error gets reduced. For the lower bound on the initial temperature, corresponding to the 
reheating temperature within inflation, the same considerations as in the SM case hold. For 
fh\ < 10 _3 eV the same lower bounds valid for M\ apply approximately also to T\. 

In the relevant strong wash-out regime for K > 4, corresponding to fh\ > 3.2 x 10 -3 eV, 
the relaxation compared to the M\ lower bound is larger than twice. In this way one obtains 
conservatively, using the upper values for Kf (cf. (|95j)). 

T { > 1.5 x 10 9 GeV. (100) 

This is an appropriate conservative value for a generical comparison with the upper bounds 
on the reheating temperature that arise by imposing that a gravitino thermal production is 
not in conflict with cosmological observations (see^J for a recent discussion and references) . 
For a more precise comparison one should calculate the lower bound on T\ for a specific 
value of rh\. For example, in the range m so i < m\ < m a t m > the relaxation compared to the 
Mi lower bound is practically the same as in the SM case (between 4 and 6 times) and thus 
one gets 

Tt > [3.5 (2) x 10 9 - 2 (1.5) x 10 10 ] GeV . (101) 

Note that these analytic results are in good agreement with the numerical ones in^4 

Let us now study how the upper bound on the neutrino masses gets modified. This can 
be easily done calculating the value of the total variation of the final asymmetry given by™ 

£ = %r^ , (102) 

where £ w is the variation of the wash-out term from AL = 2 processes that is crucial 
for the determination of the neutrino mass bound and that can be expressed through the 
parameter to in the Eq. (|59|) . We have already seen that £ e = 2. Moreover since the peak of 
the asymmetry lies in the strong wash-out regime we can also use the result Eq. 1)94(1 for £ Kf . 
Therefore, we miss only to determine £ w . There are two effects to be considered. The first 
is the increase of the number of degrees of freedom that speeds up the expansion reducing 
the efficiency of wash-out processes. The second is the presence of new different additional 
AL = 2 processes and this clearly strengthens the rate Tal=2- From the expressions given 
in 1^ one can find £r A z,=2 = 5/3 and thus £ w = £r AL=2 /\/~£g± — 5/(3 \/2). Putting all 
together one finds £ ~ 6/5 and, when the running of neutrino masses is neglected, the 



bound is about 5% more relaxed compared to the SM case, namely rnf ISSM < 0.16 eV. 
The effect of running of neutrino masses, as in the SM, goes into the direction to mak e 
the bound more stringent. However the effect can be as small as ~ 7% (for tan/3 ~ 10' 
roughly half than in the SM case, and thus one obtains in the end (at 3 a) 



mf^ M < 0.15 eV. (103) 



3.9 A 'too- short-blanket problem' 



If one requires that Mi is lower than a certain cut-off value M^ then the upper bound on 
the neutrino masses bec omes more stringent ■ * . Such a cut-off can either arise directly 
from neutrino models or indirectly from an upper bound on the reheating temperature 
T^. In this second case one has M£ ~ z*Tr where z* ~ 1 in the weak wash-out regime 
and z+ ~ zb — 2 in the strong wash-out regime k . It is then quite interesting to study 
the dependence m^ nd (M*), or equivalently m^ ou (T^). This is done in detail in "2*1 
here we just sketch some general features and results. For definiteness we will refer to the 
supersymmetric case, since in this case the avoidance of the gravitino problem implies an 
upper bound on the reheating temperature. 

First, note that if T£ < Tg 1 " 1 ~ 10 9 GeV (cf- lfTUHJl ). then simply there is no allowed 
value for mi- Moreover until one has rni° n <C m a tm ; there is a strong dependence on 
TJj, since the maximum baryon asymmetry grows linearly with T^ while is very slightly 
dependent on mi. This means that the function mj >ound (T^) has a vertical asymptote in 
TJj" n ~ 10 9 GeV. For m \ ounA ~ m a t m the suppression factor (3 ma , x (mi) = m atm /(mi + m.3) 
in the CP asymmetry and the loss in the efficiency for fh\ > m\ > 10~ 3 eV compensate the 
increase of T^ and the growth of m\ on , for increasing T^, slows down and eventually, for 
T^ > J* ca ~ 3 x 10 12 GeV, saturates to its maximum value (cf. (j!U3Jl ) and stays constant 
l . Values T^ ~ 10 10 GeV are particularly interesting since they correspond to maximum 
allowed values from gravitino problem arguments. In this case one expects mi ~ m a t m 3> 
10~ 3 eV and, since rh\ > mi, one can use the strong wash-out limit for the final efficiency 
factor Kf ~ 2/{K zb)- Moreover the wash-out factor AW from AL = 2 processes can be 
neglected. In this case it is possible to show the following approximate bound 

mi ^ A 0.2 T* . „. 

with A ~ — — —^ - , (104) 



m atm ~ VT+2A ^Catm 10 10 GeV ' 

where £„ = r/g MB /6 x lO" 10 and £ atm = m atm /0.051 eV. 

Let us consider two examples using £ v = £ atm = 1. For an upper bound T^ = 3 x 
10 10 GeV one finds m x < 0.4m atm ~ 0.02 eV, implying m 3 < 0.055 eV. For 7£ = 10 11 GeV 
one finds mi < 0.9m a t m ~ 0.045 eV, implying 7713 < 0.07 eV. 



The value of 2* has to be evaluated for that particular value of fh\ that maximizes the asymmetry. 
'Note that at the peak one has z B ^ 10 and Mf eak ~ 2 x 10 13 GeV and thus T^ ~ Mf eak /8 
3 x 10 12 GeV. 



This exercise shows that it is difficult to conciliate reheating temperatures close to the 
minimum allowed one (Tr < 10 10 GeV) and at the same time to make thermal leptogen- 
esis compatible with quasi-degenerate neutrino masses by evading the upper bound: the 
two things go into opposite directions, a typical too-short-blanket problem. There are two 
interesting consequences. The first is that if a stringent upper bound on the reheating tem- 
perature is placed, like Tr < 3 X 10 10 GeV, then it becomes difficult to evade the bound 
invoking a quasi-degenerate heavy neutrino spectrum since the bound still falls in a tran- 
sition region where light neutrinos exhibit a partial hierarchy. Vice versa if one requires 
quasi-degenerate light neutrinos to be compatible with the minimal thermal leptogenesis 
scenario then the problem of a large minimum reheating temperature gets exacerbated. 
Indeed, it is difficult in this case to avoid high values Tr > 10 11 GeV, unless one invokes 
a strong d egenera te heavy neutrino spectrum such to have a resonant enhancement of the 
asymmetry! I I Anyway further investigations are needed to understand the exact condi- 
tions on the degeneracy of the heavy neutrino spectrum and how they depend on a cut-off 
on Tr or directly on M\. 



4 Final discussion 

Leptogenesis is a specific realization of the simplest and oldest baryogenesis class of models 
where the asymmetry is generated from heavy particle decays. Its minimal version, thermal 
leptogenesis, is based crucially on neutrino properties and because of the great experimen- 
tal neutrino physics achievements it became in the last year a testable model. The decay 
parameter, the key quantity in models of baryogenesis from heavy particle decays, is a 
quantity closely related to neutrino masses. This cannot be exactly determined from data 
but there is an emerging favored range of values, K\ ep ~ 5 — 50, that implies just a small 
departure from thermal equilibrium, however large enough to explain the observed value of 
the asymmetry and moreover with some nice consequences. The predicted baryon asymme- 
try is independent on the initial conditions, both on the initial value of the asymmetry and 
on the initial number of decaying RH neutrinos. Moreover, the theoretical uncertainties are 
minimized and the final asymmetry is predicted with a precision that is within half order of 
magnitude. These features can be synthesized saying that thermal leptogenesis predictions 
are quite stable and model independent, a picture that resembles very closely the Standard 
Big Bang Nucleosynthesis in predicting the primordial nuclear abundances. The drawback 
is that values of K ~ 10 determine a loss in the efficiency between one and two orders 
of magnitude. This has to be compensated by an increase of the CP asymmetry of the 
same amount, implying a more stringent lower bound on M\. On the other hand we have 
seen how in the strong wash out regime the temperature of baryogenesis gets much smaller 
than M\ and this relaxes the lower bound on the reheating temperature compared to the 
lower bound on M\ of a factor ~ 5. Therefore, there seems to be an intriguing conspir- 
acy between neutrino mixing data and the explanation of the observed baryon asymmetry. 
Actually considerations on the maximum allowed value of the effective neutrino mass show 



that the conspiracy is even deeper^J and future experimental information on the absolute 
neutrino mass scale could give a further support. 

If the leptogenesis upper bound on the absolute neutrino mass scale, m, < 0.1 eV, will 
be fully tested with cosmology, neutrinoless double beta decay and Tritium beta decay 
experiments, then thermal leptogenesis can work in its minimal way. On the other hand, if 
neutrino masses higher than 0.1 eV will be found, then this can be either regarded as the 
effect of the existence of some level of degeneracy in the heavy neutrino spectrum, to be 
understood whether easily realized or not within the simple seesaw mechanism, or, more 
likely, as a drastic departure from the minimal thermal leptogenesis picture, at the expense 
of predictivity. In any case it should be clear, from the discussion on the 'too-short-blanket 
problem', that the two statements for which thermal leptogenesis requires dangerously large 
reheating temperatures within the supersymmetric framework and that the neutrino mass 
bound can be evaded within minimal thermal leptogenesis, can be very difficultly made 
compatible with each other. 

Another interesting aspect is that if the lightest neutrino mass m\ will be found to 
be higher than m* ~ 10~ 3 eV, then it will be possible to conclude model independently 
that thermal leptogenesis lies in the strong wash-out regime and in this way all pieces 
of the experimental information will have fitted within the theoretical best expectations. 
Therefore, if the absolute neutrino mass scale will be found to lie within the window (10~ 3 — 
10 _1 ) eV, the picture will receive further strong support from the data. 
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